The onset of the non-linear regime in unified dark matter models 
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We discuss the onset of the non-linear regime in the context of unified dark matter models involving 
a generalised Chaplygin gas. We show that the transition from dark matter-like to dark energy-like 
behaviour will never be smooth. In some regions of space the transition will never take place while 
in others it may happen sooner or later than naively expected. As a result the linear theory used in 
previous studies may break down late in the matter dominated era even on large cosmological scales. 
We study the importance of this effect showing that its magnitude depends on the exact form of the 
equation of state in the low density regime. We expect that our results will be relevant for other 
unified dark matter scenarios particularly those where the quartessence candidate is a perfect fluid. 



I. INTRODUCTION 

There is growing evidence that the Universe we live in 
is (nearly) flat and presently undergoing an accelerating 
phase 0, H H H| ■ A 'dark' energy component is thought 
to be responsible for this acceleration and either a 'stan- 
dard' cosmological constant [f| , a quintessence scalar field 
@ , or a k-essence field are usually put forth to account 
for it. However, all current explanations necessarily face 
some level of fine-tuning, which motivates a search for 
further alternatives (see for example [8ll^.lTo|). 

The Chaplygin gas provides one such interesting alter- 
native. It bears the exotic equation of state 



where p is the density, C is a positive constant and < 
a < 1. In a homogeneous universe its main property 
(cosmologically speaking) is that of mimicking cold dark 
matter at early times but progressively evolving into a 
cosmological constant later on. If a = 0, the Chaplygin 
gas model is in fact identical to a familiar ACDM scenario 
where fl^ = 1 — A is the equivalent matter density (with 
A = C/ p\ +a where p is the Chaplygin gas density at 
the present time). 

This particular behaviour indicates that dark energy 
and dark matter could just be different manifestations of 
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a unique perfect fluid known a s 'q uartessence'. In this 
regard, the Chaplygin gas EJ 0] (or simple generali- 
sations thereof EH) has attracted considerable attention 
as a quartessence prototype since a connection between 
string theory and the original Chaplygin gas has been 
proposed (see for example EH and references therein). 

The simple picture described above has been taken for 
granted in all the work on this subject so far, namely 
when confronting models with cosmological observations 
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is, however, a caveat to it which may have important 
implications. As we previously stated, the Chaplygin 
gas starts out behaving as CDM, with linear perturba- 
tions growing in the usual way as long as this regime 
holds. In linear theory, once the background evolution 
starts switching to a A-like behaviour, perturbations in 
the Chaplygin gas will become heavily damped while the 
growth of baryonic perturbations is slowed considerably. 
This would all be fine, in the sense that assuming lin- 
ear theory to hold, regions of parameter space have been 
shown to exist where structure formation proceeds in a 
way that is in agreement with observations. However, 
there is no a priori guarantee that non-linear effects will 
not significantly alter these conclusions. 

Specifically, in this letter we show that in a class of 
unified dark matter scenarios this transition from dark 
matter-like to dark energy-like behaviour will never be 
smooth. The reason is simple: at some point, density per- 
turbations on a given scale reach the non-linear regime 
and consequently a large fraction of the mass is then 
expected to be incorporated in collapsed objects whose 
evolution is effectively decoupled from the background. 
In these regions the Chaplygin gas never evolves into 
the dark energy-like stages (which only happens at suffi- 
ciently low densities). On the other hand, in low density 
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FIG. 1: The evolution of the sound speed squared, c s , of 
the Chaplygin gas as a function of redshift z for a — 0.2, 0.5 
and 1. We clearly see that at early times the Chaplygin gas 
closely resembles CDM. 
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FIG. 2: The real part of the eigenvalues of the Chaplygin 
gas. The condition Re(A(fc c ,a)) = sets the critical scale for 
that instant. The solid line represents the first eigenvalue and 
the dashed, the second. Notice the bifurcation pattern. 



regions the absolute value of the (negative) pressure may 
be much larger than the average value thus having the 
opposite effect on the dynamics of the universe (favouring 
an earlier accelerating phase). 

It is the importance of this effect that we set out to 
quantify in this letter. While the analysis will be done 
for the case of the Chaplygin gas, our results will be 
relevant for any similar scenarios where the quartessence 
candidate is a perfect fluid. We note that concerns about 
the relevance of the non-linear evolutionary phase in a 
different context (condensation and metamorphosis cos- 
mologies) have also been made in ref. |28j. 



II. THE CRITICAL SCALES 

It is straightforward to show that in a homogeneous 
and isotropic universe the Chaplygin gas energy density 
evolves as 



P = Po 



A+(l-A)(l + z) 



3(1+q) 



1/1+a 



(2) 



where A = C/pJ +a and po is the present density. In 
Fig. ^we plot the evolution of the sound speed squared 
c 2 s = aC/ p 1+a — —ap/ p as a function of the redshift z 
taking A = 0.71 for three different values of a (0.2, 0.5 
and I). At early times the sound speed squared is 



oc (1 + z) 
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(3) 



which determines how the inclination of the curves in 
Fig. n depends on a. We clearly see that deep in the 
matter era the equation of state of the Chaplygin gas 
closely resembles that of CDM (with \p/p\ -C 1). 

Using a first order perturbative analysis, the evolution 
in Fourier space, of density perturbations in the Chaply- 
gin gas, can be described by the following equation (see 



123 for details): 

5'1 + [2 + S-3(2w-c 2 s )]d' k 



(4) 



~(1 -6c 2 s +8w~3w 2 ) 



kc s 
aH 



which we rewrite as 5 k 



AS' k + BS k = 0. Here a is the 



scale factor, H is the Hubble parameter, ' = d/dlna, 
£ = H'/H and w = p/p. This second order differ- 
ential equation can easily be transformed into a non- 
autonomous dynamical system of the form x' = Ax 
where x = (5, 5') and A is the matrix of the coefficients, 
which arc functions of both the wave number k and the 
redshift z. 

We shall look for critical scales associated with the 
Chaplygin gas by calculating the eigenvalues A of A for 
fixed instants of time. These are the roots of the Cayley 
polynomial A 2 + AX + B = 0. Since the sign of the real 
part of the eigenvalue determines whether or not a mode 
grows or decays, the condition Re(A(fc c ,a)) = sets the 
critical scale for that a. Fig. 2 contains the real part of 
the two eigenvalues for several times. A bifurcation pat- 
tern emerges. At early times the Chaplygin gas acts as 
CDM and so all relevant scales are gravitationally unsta- 
ble. As time goes by, this critical scale starts to diminish 
until there are no gravitational unstable scales. 

We note that in this linear analysis both c s and w are 
evaluated to zeroth order. We will show that this is no 
longer a good approximation when the fluctuations in the 
Chaplygin gas component become large. 



III. NON-LINEAR EFFECTS 

Having gained some intuition for the relevant length 
and time scales in the problem, we now proceed with 
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FIG. 3: The linear evolution of the dispersion, a, of the den- 
sity fluctuations in the baryon (solid line) and Chaplygin gas 
(dashed lines) components as a function of R and a assuming 
a = 1. Note that at early times the baryon and Chaplygin 
gas fluctuations evolve in the same manner. Later on, pres- 
sure effects prevent the Chaplygin gas from collapsing further. 
However the baryon fluctuations can still keep growing (at a 
slower pace). 



a more detailed analysis. Throughout the analysis we 
adopt priors in agreement with the WMAP first-year 
data release with a power-law initial spectrum 4]. The 
parameters are an equivalent matter density fl^ = 1 — 
A + fib = 0.29, a baryon density Qb = 0.047, an equiva- 



lent cosmological constant density f2^ 



1-0! 



0.71, 



a Hubble parameter h — 0.72, a normalisation ag — 0.9, 
and perturbation spectral index n s = 0.99. 

We start by determining how the linear density per- 
turbations of the Chaplygin gas and baryon components 
evolve with redshift (the details of such an analysis are 
described in Sect. II of our earlier work (2(J). We empha- 
sise that the use of linear theory early in the matter era 
is actually a good approximation on large enough scales. 
The effects of the breakdown of linear theory will only 
be important (on large cosmological scales) at late times 
when a smooth transition from a dark matter to a dark 
energy dominated universe would naively be expected to 
take place. 

Hence, we use linear theory in order to compute the 
value of the dispersion of the density fluctuations in the 
baryon and Chaplygin gas components, a(R,a), as a 
function of R and a. This is plotted in Fig. [3] for the 
particular case of the original Chaplygin gas with a — 1 . 
We see that since the Chaplygin gas behaves as matter 
at early stages (see Fig. QJ, embedded perturbations will 
grow proportionally to the scale factor, a, and in tune 
with those in the baryonic component. Therefore, both 
fluids evolve in the same way earlier on and have ap- 
proximately the same value of a on all relevant scales. 
During this stage we have a oc a — (1 + z)^ 1 . Later 
on, the pressure of the Chaplygin gas will have increased 



FIG. 4: The ratio between the average value of p and and 
its zeroth order background value, (p)/pb, as a function of 
the dispersion, a, of the linear density fluctuations in the 
Chaplygin gas component for various values of a (1, 0.8, 0.6, 
0.4, 0.2 and from top to bottom). For a > we clearly see 
that {p)/pb rapidly diverges from unity if a is large enough. 



dramatically preventing further collapse of the Chaply- 
gin gas component. However, the baryon fluctuations 
can still keep growing (at a slower pace). We also see 
in Fig. [3] that the Chaplygin gas component becomes 
non-linear on small scales early in the matter era. It is 
clear that when this happens a significant fraction of the 
Chaplygin gas will have collapsed and decoupled from the 
background so that a transition from a dark matter-like 
to a dark energy-like stage (which necessarily requires 
lower densities) never happens in those regions. Using 
the Press- Schechter framework [2^| we can show that for 
cr = l the fraction of the equivalent mass that is already 
incorporated in collapsed objects is already close to 0.1. 

In the non-linear regime an initial gaussian density 
field is better described by a lognormal one-point prob- 
ability distribution function, V{8), given by (see for ex- 
ample [30j and references therein) 



V(S) 



(1 + 5)- 1 



exp 



in 2 ((i+^yr+^y 

21n(l + ^) 



(5) 

where a 2 ^ = exp(cr 2 ) — 1 and a is computed using linear 
theory. We use 10) in order to calculate the ratio between 
the average value of p and its zeroth order background 
value (p b = -C/p%), 



(p)/Pb 



(l + S)- a V(S)dS. 



(6) 



as a function of the dispersion of the density fluctua- 
tions in the Chaplygin gas component, a. The results 
are displayed in Fig. 0] for various values of a (0, 0.2, 
0.4, 0.6, 0.8 and 1). We see that in all but one case 
(for a = 0) the average value of the pressure strongly 
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diverges from its zeroth order background value as soon 
as a becomes large enough. This means that the average 
values of both w and c 2 s will also move away from their 
zeroth order value (except in the a = case) causing the 
breakdown of linear theory. The magnitude of this effect 
becomes more pronounced at late times when the nega- 
tive pressure starts to become dynamically important on 
all scales. 

We further note that in the a / case and for very 
small densities the sound speed, c s oc / o _ ( 1 + Q )/ 2 j will be- 
come greater than the speed of light, c = 1. Hence, we 
may expect that in the context of realistic models there 
will be a cut-off to this power law behaviour. One equa- 
tion of state that resembles that of eqn. at high den- 
sities while switching off to a cosmological constant like 
equation of state at low densities is 



P = 



P 



(P + A) 



1+/3 



(7) 



where A > is a low density cut-off and Q < j3 < 1. Note 
that in this case the absolute values of both c 2 s and w are 
always smaller than unity. Another interesting case is 



P 



C 



(p + A)' 3 



(8) 



for which c 2 s < 1 (assuming that C/3/A /3+1 < 1) but the 
absolute value of w can grow arbitrarily large. As we have 



shown the importance of the non-linear effects discussed 
in this letter for cosmological predictions will strongly 
depend on the exact form of the equation of state in the 
low density regime. 



IV. CONCLUSIONS 

We have shown that previous work on the Chaplygin 
gas has a caveat which may have crucial implications for 
the predicted observational consequences of the model. 
The linear theory used in previous treatments breaks 
down at late times even on large cosmological scales ex- 
cept in the a — case. This means that non-linear effects 
should be taken into account when confronting the model 
with cosmological observations. 

Although the present analysis has been done for the 
particular case of a Chaplygin gas we expect our results 
to be relevant for other unified dark matter scenarios 
particularly those where the quartessence candidate is a 
perfect fluid. 
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